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4.  Development  of  the  Novel  Simulation  Approach 


This  project  developed  a  new  simulation  scheme  for  investigating  electromagnetic  wave 
effects  on  millimeter-wave  transistors  and  other  semiconductor  devices.  This  scheme  is  based  on 
coupling  a  set  of  equations  capable  of  accurately  describing  the  electron  transport  phenomena  in 
semiconductors,  with  a  three-dimensional  solution  of  Maxwell’s  equations.  At  this  stage,  the 
active  device  model  is  based  on  the  moments  of  the  Boltzmann's  transport  equation  obtained  by 
integration  over  the  momentum  space.  Three  equations  are  used;  Carrier  conservation  equation, 
Momentum  Conservation  equation,  and  Energy  conservation  equation.  These  equations  are 
coupled  with  a  full-electromagnetic  wave  model.  The  complete  set  of  equations  are  solved  in 
time  domain  after  discretizing  it  over  a  three  dimensional  mesh  extending  over  the 
semiconductor  material,  the  surrounding  media,  and  the  electrodes  including  the  distribution 
pads.  The  coupling  between  the  two  models  is  established  by  using  fields  obtained  from  the 
solution  of  Maxwell’s  equations  in  the  active  device  model  to  calculate  the  current  densities 
inside  the  device.  These  current  densities  are  used  to  update  the  electric  and  magnetic  fields. 

The  funding  for  this  program  started  in  May,  1995.  The  tasks  accomplished  in  the  first 
year  can  be  summarized  as  follows: 

-  The  code  for  the  obtaining  the  semiconductor  device  characteristics  was  completed  and 
adjusted  for  coupling  to  the  full-wave  solver. 

-  The  code  for  the  Electromagnetic  solution  was  developed  and  adapted  to  accept  the 
semiconductor  device  equations. 

-  The  two  models  were  successfully  merged. 

-  Elaborate  system  of  absorbing  boundary  conditions,  nonuniform  mesh,  electrode  pads  and 
various  other  improvements  were  introduced. 

The  complete  code  was  used  to  investigate  several  electromagnetic  waves  and 
semiconductor  devices  interactions. 

5.  Significance  of  the  New  Simulator  Development 

We  completed  the  development  of  the  FIRST  high  frequency  device  simulator  that 
combines  the  electromagnetic-waves  with  the  electron  transport  for  .  This  is  a  physically-based 
model,  which  capable  of  representing  short  gate  effects  and  other  aspects  that  make 
semiconductor  devices  operate  at  high  speed/high  frequencies,  as  well  as  the  effects  related  to 
electromagnetic  wave  propagation.  It  will  be  referred  to  as  a  combined  electromagnetic  and 
solid-state  (CESS)  simulator.  Usually,  microwave  devices  are  modeled  as  active  devices 
embedded  between  two  ideal,  lossless  transmission  lines  or  in  relatively  simple  circuit  models. 
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Such  approaches  has  two  main  drawbacks:  1)  Quasi-TEM  transmission  lines  are  used,  and  2)  the 
semiconductor  device  model  is  developed  in  total  separation  from  the  EM  wave.  Therefore, 
direct  solutions  of  Maxwell's  equations  are  needed  for  a  more  accurate  and  general  approach. 

The  full-wave  physical  model  used  in  this  work  allows  a  flexible  description  of  the  device  along 
with  the  appropriate  representation  of  simulation  parameters.  The  application  of  the  model  to 
different  electronic  structures,  such  as  the  High  Electron  Mobility  Transistor  (HEMT),  with 
different  material  profiles  and  boundaries  is  straight  forward.  This  was  accomplished  and 
presented  in  our  publications  ( please  see  sample  of  publications). 

6.  Simulation  of  a  Complete  Amplifier 

After  developing  the  CESS,  our  next  goal  was  to  simulate  a  microwave  amplifier  with  the 
passive  matching  networks  using  the  CESS  model.  The  simulation  of  a  complete  amplifier  using  the 
Finite-Difference  Time-Domain  (FDTD)  method  requires  intensive  computer  memory  and  consumes 
a  considerable  amount  of  time.  The  problem  of  large  computer  memory  and  time  can  be  reduced  by 
breaking  the  amplifier  circuit  into  active  and  passive  parts,  then  model  the  passive  parts  separately 
using  the  FDTD  technique  and  combine  them  with  the  full-wave  simulation  of  the  transistor. 
Breaking  a  large  circuit  into  sub-circuits  can  be  achieved  by  using  the  time-domain  diakoptics 
method  or  the  discrete  time-domain  Green’s  function. 

The  microwave  amplifier  consists  of  input  and  output  matching  networks  and  a  transistor.  The 
input  and  the  output  matching  networks  are  very  large  compared  to  the  transistor.  But  the  mesh  size 
and,  consequently,  the  FDTD  stability  criteria  for  the  amplifier  are  severely  limited  by  the  debye 
length  of  the  semiconductor.  This  imposes  a  constraint  on  the  time  step  At  of  the  FDTD  algorithm, 

which  becomes  on  the  order  of  10' 17  seconds.  If  the  input  and  the  output  matching  networks  were  to 
be  simulated  using  this  criteria,  it  would  take  so  long  time  to  simulate  the  amplifier  circuit  that  it 
becomes  almost  impossible  using  the  current  computers.  In  this  research,  a  fiill-wave  analysis  is 
performed  to  simulate  the  microwave  amplifier  with  two  tuned  coplanar  wave  guides  (CPW)  as 
matching  networks  (see  Fig.  1).  The  whole  amplifier  is  divided  into  three  regions.  The  simulation  of 
each  region  is  performed  individually  and  coupled  to  the  next  stage  properly  with  all  the  required 
information  from  the  preceding  stage.  This  technique  enables  one  to  use  large  space  step,  and  hence, 
large  time  step  in  matching  networks.  The  computer  simulation  time  is  reduced  drastically  compared 
to  the  method  incorporating  the  non-uniform  mesh  for  the  whole  amplifier.  The  computer  memory 
requirement  is  also  lowered  by  approximately  66  %  at  a  certain  time. 

Our  group  accomplished  this  major  milestone  in  the  process  of  developing  physical 
models  for  microwave  circuits  in  1997.  It  should  be  noted  that  this  is  the  first  time  a  complete 
amplifier  is  modeled  using  the  physical-based  simulations.  This  is  the  only  simulation  relating 
the  high  frequency  amplifier  performance  to  the  microscopic  electronic  characteristics  of  the 
active  device  and  the  electromagnetic  matching  network  simultaneously.  Several  research 
groups,  in  US  and  abroad,  are  developing,  or  attempting  to  develop,  similar  simulators  following 
our  approach. 
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Fig.  1.  GaAs  transistor  amplifier  with  CPW  matching  networks  (a)  Schematic  (b)  Detailed 
diagram. 


7.  Comprehensive  Analysis  of  Discontinuity  Effects  on  High  Frequency 
Transistors 

In  conventional  high  frequency  GaAs  transistors,  several  configurations  are 
implemented.  The  layout  of  the  physical  device  affects  its  performance,  especially  at  high 
frequency  of  operation.  Typical  configurations  for  a  MESFET  or  HMET  are  the  interdigitated 
construction  and  the  four-finger.  It  is  clear  that  the  electrodes  of  the  active  device  form  a  guided 
wave  structure.  At  high  frequency,  those  electrodes  have  to  be  treated  as  transmission  lines.  This 
transmission  lines  can  be  regarded  as  slot  lines  or  coplanar  waveguide,  which  forms  a 
discontinuity  at  the  input  and  output  ports  of  the  device.  Therefore,  the  input  and  output 
impedance  must  be  carefully  analyzed.  A  common  practice  in  analyzing  the  discontinuity  effects 
on  active  devices  is  to  analyze  the  passive  structure  separately  from  the  behavior  of  the  real 
active  device.  This  approach  does  not  take  into  account  the  interaction  between  the  wave 
propagating  along  the  device  electrodes  and  the  electron  transport  in  the  active  device  and  then- 
effects  on  the  device  characteristics.  On  the  other  hand,  it  is  very  difficult  to  determine  the 
discontinuity  effects  from  measurement.  This  situation  prevents,  in  many  cases,  the  possibility  of 
matching  the  FET  at  all.  Consequently,  the  true  performance  of  the  FET  might  be  unknown. 
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In  this  research,  an  algorithm,  that  deals  with  the  electromagnetic  wave  and  the 
electrodynamics  inside  the  transistor  simultaneously,  is  used  to  study  the  discontinuity  effects  on 
the  high  frequency  performance  of  MESFET  transistor.  Although  the  method  is  difficult  to 
formulate,  it  allows  more  accurate,  flexible  and  realistic  conditions  to  be  investigated.  The 
physical  semiconductor  model  is  based  on  the  complete  hydrodynamic  system.  The  full-wave 
solution  is  used  to  update  the  ac  fields  inside  the  device.  These  two  schemes  are  blended  together 
using  the  current/field  relation.  The  s-parameters  are  calculated  for  different  device 
configurations.  Very  informative  results  regarding  the  effects  of  discontinuities  on  millimeter- 
wave  transistor  performance  were  generated.  They  are  useful  for  device  optimization  and 
developing  new  designs  that  exploits  the  wave  propagation  effects. 

8.  Extending  the  applications  to  quasi-optical  power  combiners 

The  power  and  the  versatility  of  the  FDTD  technique  has  been  limited  by  its 
computational  cost,  specially  when  it  comes  to  the  simulation  and  optimization  of  circuits  with 
large  areas.  In  addition,  the  method  has  accuracy  and  stability  problems  when  simulating 
structures  which  require  variable  mesh  coarseness,  like  in  the  case  of  simulating  circuits  with 
active  devices.  The  diakoptics  technique  overcomes  these  difficulties  by  breaking  down  large 
structures  into  smaller  ones  and  characterizing  each  subsection  with  its  impulse  response.  The 
analysis  of  the  full  structure  is  then  achieved  by  subsequently  interconnecting  the  individual 
responses  together  via  a  convolution  process.  This  method  is  well  established  for  the  analysis  of 
steady  state-networks. 

The  purpose  of  this  research  is  to  expend  and  validate  time  domain  diakoptics  to  3D  structures  to 
extend  this  approach  to  millimeter-wave  quasi-optical  systems.  Techniques  to  save  memory  and 
computational  time  are  investigated  and  developed.  Validation  of  this  technique  is  illustrated  by 
comparing  the  results  obtained  using  this  method  with  those  obtained  by  using  FDTD  simulation 
of  the  full  structure.  Excellent  agreement  was  obtained.  This  technique  was  applied  to  the 
antenna  used  in  a  two-dimensional  quasi-optical  power  combiner.  The  objective  is  to  model  the 
entire  structure,  including  the  active  devices,  using  the  combined  electromagnetic  solid-sate 
simulator. 

9.  Development  of  the  3-D  diakoptic  Approach  for  circuit  simulations 

The  power  and  the  versatility  of  the  FDTD  technique  has  been  limited  by  its  computational  cost, 
specially  when  it  comes  to  the  simulation  and  optimization  of  circuits  with  large  areas.  In 
addition,  the  method  has  accuracy  and  stability  problems  when  simulating  structures  requiring 
variable  mesh  coarseness,  like  in  the  case  of  simulating  circuits  with  active  devices.  The 
diakoptics  technique  overcomes  these  difficulties  by  breaking  down  large  structures  into  smaller 
ones  and  characterizing  each  subsection  with  its  impulse  response.  The  analysis  of  the  full 
structure  is  then  achieved  by  subsequently  interconnecting  the  individual  responses  together  via  a 
convolution  process.  Our  contribution  to  this  subject  is  two  fold:  First,  we  expended  and 
validated  time  domain  diakoptics  for  3-D  microwave  structures  (validation  of  this  technique  is 
illustrated  by  comparing  the  results  obtained  using  this  method  with  those  obtained  by  using 
FDTD  simulation  of  the  full  structure).  Second,  We  used  the  FDTD/diakoptics  method  to 
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perform  a  full-wave  simulation  of  the  passive  structure  of  quasi-optical  amplifier  used  in  two- 
dimensional  power  combining  systems.  To  limit  the  memory  requirement  the  convolution 
process  is  performed  only  at  few  selected  points  of  the  port  mesh.  The  field  distribution  at  the 
rest  of  the  mesh  points  is  approximated  by  modulating  the  calculated  fields  with  the  steady  state 
field  distribution.  To  reduce  the  heavy  burden  of  the  convolution  operation,  the  Prony  method  is 
used  to  approximate  the  impulse  responses  and  generate  a  fast  recursive  convolution.  Details  of 
this  method  can  be  found  in  our  published  papers  listed  hereafter. 

10.  Low  Loss  Air-Gap  Transmission  Lines  and  Spiral  Inductors  for  MMICs  Using  Glass 
Microbump  Bonding  Technique 

One  of  the  byproducts  of  our  research  is  to  use  knowledge  and  insight  gained  to  develop  better 
active  and  passive  devices  for  high  frequency  operations.  For  example,  new  types  of 
transmission  lines,  called  air-gap  lines,  were  developed,  simulated,  fabricated  and  tested  in  our 
group.  These  are  microwave  transmission  lines  deposited  on  glass  superstrates  and  then 
connected  to  active  devices,  fabricated  on  lossy  substrates  as  usual,  using  a  new  technology 
known  as  Glass  Microbump  Bonding  Technique  GMBB  (also  developed  in  our  group.)  The  low- 
loss  characteristics  is  one  of  the  main  advantages  of  the  new  air-gap  transmission  lines. 
Moreover,  air-gap  spiral  inductor  structures  have  been  fabricated  and  integrated  with 
semiconductor  substrates  using  the  same  technology.  Spiral  inductors  using  air-gap  structures 
have  the  advantages  of  low  losses,  and  low  parasitic  capacitance  compared  to  conventional 
inductors  on  doped  silicon  semiconductor  substrate.  Stacked  air-gap  spiral  inductors  on  GaAs 
substrates  using  GMBB  techniques  also  can  reduce  the  inductor  area.  Because  the  glass 
microbump  bonding  techniques  are  simple,  this  bonding  technique  provides  an  alternative 
integration  approach  for  monolithic  microwave  integrated  circuits  (MMICs).  Experimental 
results  of  air-gap  spiral  inductor  on  both  silicon  and  GaAs  substrates  were  obtained. 

11.  Technology  Transfer 

The  interest  in  this  novel  technique  for  analyzing  semiconductor  devices  prompted  us  to 
widely  disseminate  the  information.  More  than  50  archival  journal  papers,  dissertations,  and 
presentations  has  resulted  from  this  research.  There  are  listed  in  the  following  sections. 


All  the  personnel  funded  from  this  project  has  joined,  or  accepted  offer,  from  industrial 
companies  in  US  ( The  P.I.  is  not  included).  Their  names  and  current  employers  are  as  follows: 

1.  Mohamed  A.  Megahed,  Ph.D.,  Post  Doctoral  Research  Associate,  joined  “Peregrine 
Semiconductors.” 

2.  S.M.  Sohel  Imtiaz,  Ph.D.  Graduate  Research  Assistant,  joined  “Micro  Linear.” 

3.  JefFCheuang,  Ph.D.,  Graduate  Research  Assistant,  joined  “M/A-COM,  Inc.” 

4.  Samir  Hammadi,  Ph.D  expected  in  May  99.  Graduate  Research  Assistant,  Has  already 
received  some  offers  from  companies  in  US,  decision  pending. 
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(a)  Side  View 


lass  Microbump 


(b)  Top  View 


Fig.  2.  Air-gap  microstrip  line  (AGML)  structure  fabricated  using  GMBB  techniques. 


(a) 


(b) 


(c) 


Fig.  3.  Photographs  of  various  types  of  inductors,  (a)  Conventional  inductor,  (b)  Air-gap 
inductor,  (c)  Stacked  inductor. 


12.  Conference  and  Workshop  Presentations 

[1]  M.  Alsunaidi  and  S.  El-Ghazaly,  “Electromagnetic  Wave  Effects  on  Millimeter  Wave 
Semiconductor  Devices;”  PIERS  19952  Progress  in  Electromagnetic  Research 
Symposium,  Seattle,  WA ,  P.  795,  July  1995. 

[2]  R.  Grondin,  I.-S.  Lim  and  S.  El-Ghazaly,  “Parallel  Processing  of  Full-Wave 
Electromagnetic  Models  of  Ultrafast  Photoconductive  Switching;”  PIERS  1995- 
Progress  in  Electromagnetic  Research  Symposium,  Seattle,  WA,  P.  801,  July  1995. 

[3]  S.M.  Sohel  Imtiaz,  M.A.  Alsunaidi,  M.A.  Megahed  and  S.  El-Ghazaly,  "Electromagnetic 
Wave  Effects  on  MODFETs;"  1996  National  Radio  Science  Meeting -  U.S.  National 
Committee  ofURSI,  Boulder,  Co,  p.  155,  January  1996. 
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16.  Special  Effort  in  Information  Dissemination 

Besides  this  long  list  of  conference  and  journal  papers,  the  P.I.  devoted  a  special  attention 
to  widely  disseminate  the  information  and  the  results  generated  from  this  research.  This  is  a 
particularly  important  aspect  of  this  project  given  its  novelty,  which  necessitates  that  the 
wide  majority  of  the  microwave  and  RF  community  need  to  be  informed  about  it.  This  effort 
included  organizing  special  sessions  in  international  conferences,  inviting  speakers, 
organizing  panel  discussions,  and  workshops.  This  wide  dissemination  effort  resulted  in  a 
special  invitation  to  become  an  editor  of  a  special  issue  of  the  IEEE  Transactions  on 
Microwave  Theory  and  Techniques  on  this  topic,  with  two  other  colleagues.  This  special 
issue  will  appear  in  June  1 999. 

It  has  been  observed  that  several  groups  all  over  the  world  recognized  the  importance  of 
this  approach  and  directed  their  research  to  it.  Continuously  increasing  number  of  conference 
and  journal  papers  are  being  published.  It  is  only  a  matter  of  time. 
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Amplifiers 

3.  Integration  of  Air-Gap  Transmission  Lines  on  Doped  Silicon  Substrate  Using  Glass 
Microbump  Bonding  Techniques 

4.  Parallel  Implementation  of  a  Two-Dimensional  Hydrodynamic  Model  for  Microwave 
Semiconductor  Devices  Including  Inertia  Effects  in  Momentum  Relaxation 
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